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ABSTRACT 

We report the first results from evolutions of generic black-hole binaries, i.e. a binary containing 
unequal mass black holes with misaligned spins. Our configuration, which has a mass ratio of 2 : 1, 
consists of an initially non-spinning hole orbiting a larger, rapidly spinning hole (specific spin a/m = 
0.885), with the spin direction oriented —45° with respect to the orbital plane. We track the inspiral 
and merger for ~ 2 orbits and find that the remnant receives a substantial kick of 454kms~ , more 
than twice as large as the maximum kick from non-spinning binaries. The remnant spin direction is 
flipped by 103° with respect to the initial spin direction of the larger hole. We performed a second run 
with anti-aligned spins, a/m = ±0.5 lying in the orbital plane that produces a kick of ~ 1830kms _1 
off the orbital plane. This value scales to nearly 4000 km s~ for maximally spinning holes. Such a 
large recoil velocity opens the possibility that a merged binary can be ejected even from the nucleus 
of a massive host galaxy. 

Subject headings: black hole physics - galaxies: nuclei - gravitational waves - relativity - gravitation 



1. INTRODUCTION 

One of the major goals of numerical relativity has 
been the accurate evolution of generic black-hole bina- 
ries from inspiral through merger and ringdown. It is in 
this non-linear merger regime where most of the gravita- 
tional radiation is emitted, including the radiation of lin- 
ear momentum responsible for large merger recoils that 
can eject the remnant from the host galaxy. With the 
recen t breakthroughs in numer i cal techniques (iPretoriusl 
120051: ICampanelli et al.l l2006at iBaker et all l2006aT ) this 
goal is finally being realized. Within the past 18 
months rapid progress has been achieved in our un- 
derstanding of black-hole-binary mergers. Non-spinning 
equal-mass binarie s were studied in deta i l, including 
the las t few orbits (|Campanelli et al. 2006b; Bakc r~et al.l 
I2006U) . the effects of elliptical motion (|Pretoriusll2006f) 
on the gravitational radiation, and waveforms gener- 
ated from binaries with large initial separations were 
successfully matched to post- Newtonian theory with 
very goo d agreement (jBuonanno et al . 2006; Baker et al.l 
Non-spinnin g uneq ual mass black holes were 



studied in ICampanellil (120051). I Herrmann et all (|2006f) . 
IBaker et al.l (|2006eD and IGonzalez et al.l (|2006T ). where 
the recoil velocity of the post-merger remnant was 
co mputed. In particu lar, the accurate calculations 
of IGonzalez et al. (2006) indicate that the maximum re- 
coil velocity of non-spinning quasi-circular binaries, with 
mass ratio q — m\/mi w 1/3, is ~ 175 kms 1 . Simu- 
lations of highly-spinning b l ack- ho le binaries were intro- 
duced in ICampanelli et all (|2006ch where it was shown 
that the direction of the spin (in that case either aligned 
or counter-aligned with the orbital angular momentum) 
has a strong effect on the merger time and energy- 
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momen tum radiated to infinity. In Campanell i et all 
(2006d) it was found that the non-linear tidal effects 
were too weak to drive a binary int o a corotating state. 
Finally in ICampanelli et alJ (|2006el ) spin precession and 
spin-flips were studied for equal-mass binaries with indi- 
vidual spins not aligned with the orbital angular momen- 
tum (but with individual spins having the same magni- 
tude and direction). 

All the previous simulations contained symmetries 
which suppressed some important astrophysical proper- 
ties (e.g. precession, recoil, spin-orbit coupling) of generic 
binary mergers, and in the case of the recoil calculation, 
spins were neglected entirely. With the knowledge gained 
from these simulations we can now design and evolve a 
truly prototypical black- hole binary. In the scenario con- 
sidered here, a high-mass black hole, with a specific spin 
of a/m — 0.885 (the largest considered thus far), merges 
with a smaller hole having negligible spin. The mass ratio 
of the two holes is 1.99, and the initial binary configu- 
ration is such that the spin of the larger hole points 45° 
below the orbital plane. This configuration will manifest 
precession of the spin axis, a significant spin flip of the 
remnant spin with respect to the initial individual hori- 
zon spin, and a significant recoil kick. The simulations 
that we report in this letter show that the recoil due to 
the spin can be more than an order of magnitude larger 
than the maximum recoil due to unequal-masses alone. 

2. TECHNIQUES 

We use the puncture approach ([Brandt and Brugmann 
1997 ) along with the TwoPunctures ( [Ansorg et al. 
20041 ) thorn to compute initial data. We evolve 
these b lack hole binary data sets using the 
LazEv ([Zlochower et al.l I2005D implementation of the 
moving puncture approach (ICampanelli et all l2006aT) 
(which is based on th e BSS N dNakamura et all 119871: 
Shibata and Nakamural 119951 : iBaumgarte and Shapiro ! 
1999) formulation). We use the Carpet (|Schnetter et al.l 
20041 ) mesh refinement driver to provide a 'moving 
boxes' style mesh refinement. In this approach refined 
grids of fixed size are arranged about the coordinate 
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TABLE 1 

Initial data parameters for the SP6 (top) and SP2 

(BOTTOM) CONFIGURATIONS. m p IS THE PUNCTURE MASS 
PARAMETER OF THE TWO HOLES. SP6 HAS SPINS Si = (0, S, —S) 

AND S2 = (0,0,0), MOMENTA P = ±(P r , P ± ,0), PUNCTURE 
POSITIONS XI = (x+,d,d) AND 22 = (x-,d,d), AND MASSES mi 

and m.2. While SP2 has spins Si = — S2 = (0, S, 0), puncture 

POSITIONS Xi = —X2 = (a:, 0,0), AND momenta 
Pi = -P 2 = (0,P,0). 



mp/M 


0.3185 


d/M 


0.0012817 


mi/M 


0.6680 


x + /M 


2.68773 


P r /M 


-0.0013947 


m 2 /M 


0.3355 


x-/M 


-5.20295 


P±/M 


0.10695 


S/M 2 


0.27941 


■trip/M 


0.430213 


x/M 


3.28413 


S/M 2 


0.12871 


P/M 


0.13355 


m/M 


0.5066 







centers of both holes. The Carpet code then moves these 
fine grids about the computational domain by following 
the trajectories of the two black holes. We measure 
the horizon spin (ma gnitude and direction) u sing the 
techniques detailed in ICampanelli et al.l (|2006ef h 

3. RESULTS 

The initial data parameters for our SP6 configuration 
(i.e. generic binary configuration), which were obtained 
using the 3PN equations of motion, are given in Table [TJ 
Note that the binary has a small inward radial velocity 
(that we obtain from the post-Newtonian inspiral.) The 
initial orbital plane coincides with the xy plane. 

We tested our code wi th mesh refinemen t by ev olving 
the SP3 configuration of ICampanelli et all |2P06c). For 
this test we evolved SP3 with three different grid con- 
figuration with finest resolutions of M/32, M/40, M/52 
respectively, and 6 levels of refinement. We placed the 
refinement boundaries at the same coordinate distance 
from the punctures for each configuration. We confirmed 
that the waveforms converge to fourth order and agree 
with our unigrid SP3 evolution. 

We ran the SP6 run with 7 levels of refinement, with a 
finest resolution of M/43.6. The outer boundaries were 
placed at 250M. We tracked the individual horizon spins 
throughout the evolution and found no significant spin- 
up of the smaller (initially non-spinning) hole (the value 
of a/m at merger was ~ 10~ 4 ). The larger black hole, on 
the other hand, showed a significant 45° angle of spin- 
precession, with final spin (at merger) S™ erger /M 2 = 
(-0.262, 0.189, -0.214). During the evolution the binary 
performed ~ 1.8 orbits prior to the formation of the com- 
mon apparent horizon (CAH). The first common horizon 
was detected at t C An/M = 197.96 ± 0.07. The common 
horizon had mass M H /M = 0.9781 ± 0.0001 indicating 
that (2.19 ± 0.01)% of the mass was converted into grav- 
itational radiation. The spin of the remnant horizon was 
S Iem /M 2 = (-0.0397 ± 0.0005,0.242 ± 0.002,0.4097 ± 
0.0002). The ADM angular momentum for this system is 

Jadm/M 2 = (0,0.27941,0.56447), thus we predict that 
the radiated angular momentum is J ra d/M 2 = (0.0397 ± 
0.0005, 0.037 ±0.002, 0.1548 ±0.0002). The measured ra- 
diated mass and angular momentum, based on the t — 2 
through t = 4 modes of ip 4 were E rad /M = 0.0218 ± 

0.0004 and J rad = (0.04 ± 0.01, 0.04 ± 0.01, 0.16 ± 0.01), 
which agree well with the remnant horizon parameters. 
Note that the agreement in J rac i between the horizon 
spin and radiation calculation indicates that our method 




Fig. 1. — The recoil velocities for the SP6 configurations as 
measured for an observed at r = 30 M . 



for calculating the spin direction is reasonably accurate 
for our choice of coordinates. The final horizon spin is 
flipped by 103° with respect to the initial spin of the 
larger individual horizon and 33° with respect to the ini- 
tial orbital angular momentum. 

The remnant hole acquires a significant recoil velocity 
of vLoii = (-208 ± 30, -48 ± 7,424 ± 10)kms~ 1 (see 
Fig-QJj which makes an angle of 27° with respect to the 
initial orbital angular momentum and 135° with respect 
to the initial spin. We measured this kick by calculating 
the r adiated linear momentum (Campan elli and Loustol 
1999) based on the 1 — 2 through i = 4 modes of ^4- 
We extracted these modes at r = 25M, 30M, 35M, 40M 
and then extrapolated the radiated momenta calculated 
at these radii to r = 00 using a linear (least-squares) fit 
(we excluded the initial data burst from this momentum 

calculation). The quoted errors in Vrecoii are the differ- 
ences between the linear extrapolation and a quadratic 
extrapolation. This recoil velocity of 454 ± 25kms~ 1 is 
more than double t he maximum recoil ve locity found for 
non-spinning holes (Gonzal ez et al.ll2006l) even including 
small eccentricity effects (jSopuerta et al.ll2006t) . Further- 
more, the spin-induced recoil in the xy plane might be 
offset by the mass-diffcrencc-induced recoil, potentially 
implying that a rotation of the spin about the z-axis may 
lead to a significantly larger in-plane component of the 
recoil velocity. Further study will be needed to determine 
if the mass-difference-induced recoil is (partially) aligned 
or counter-aligned with the spin-induced recoil. 

4. DISCUSSION 

We studied, for the first time using fully non-linear nu- 
merical relativity, a realistic astrophysical configuration 
of unequal mass, spinning black holes starting from a 
slightly elliptical orbit, with radial inward velocity as pre- 
dicted by post-Newtonian theory, for large initial separa- 
tions. Our main new result is that the spin component to 
the recoil velocity may produce the leading contribution. 
This is suggested by the fact that the z-component of 
the recoil velocity, which is not present for non-spinning 
binaries, is the dominant component. This also can be 
seen from the 2nd post-Ne wtonian expr essions for the 
radiated linear momentum (|Kidderlll995ti 

'-■ 8 /i 2 m 



. _ 5 : |4r(vx A) - 2v 2 (h x A) 
-(hxv) 3r(n ■ A) + 2(v ■ A) j, 



(1) 
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where x = x\ — X2, r = |x|, v = dx/dt, h = x/r, (J, = 
m\mil"m, m = mi + rri2, A = ra{Sijr(ii — Si /mi), and 
an overdot denotes d/dt. 

Based on this expression we can predict that the maxi- 
mum recoil velocity is reached for equal mass black holes 
with opposite (and maximal) spins lying on the orbital 
plane since all four terms add constructively to the ra- 
diated momentum. We performed one additional run, 
denoted by SP2, with anti-aligned spins of magnitude 
a/m ± 0.5 lying initially along the y-axis as reported in 

Table U We obtain a V lecoii = (0,0, 1830 ± 30)kms _1 . 
By rescaling this to maximal spins we obtain essentially 
double those values raising the maximum recoil of spin- 
ning holes to almost 4000 km s -1 . 

Equation (JXJ) also allows us to propose an empirical 
formula for the total recoil velocities 4 

V IeC oil(q,ct) = v m ei +uj_(cos£ei + sin£e 2 ) +v\\ e z , 



5. ASTROPHYSICAL IMPLICATIONS 



(i + <z) s 



V± = H JTTW ~' y "''- 



v {l =Kco S (e - 9 ) (i ^ g)5 (ai - qai) , (2) 

where on = Si/mf, the index _L and || refer to perpen- 
dicular and parallel to the orbital angular momentum re- 
spectively, ei, e% are orthogonal unit vectors in the orbital 
plane and £ measures the angle between the unequal mass 
and spin contribution to the recoil velocity in the orbital 
plane. The constants H and K can be determined from 
newly available runs. The angle is defined as the angle 
between the in-plane component of A and the infall direc- 
tion at merger. We have confirmed this cos O dependence 
by evolving a set of runs similar to SP2 but with initial 
spins rotated by 5<d = 7r/4, it/2, and tt. The resulting 
kicks were well modeled by a cos(0 — 0o) dependence, 
with 14 = (1873±30)kms- 1 cos[^e-(0.18±0.02)]. Note 
that we measured a maximum kick of (1830 ± 30)kms _1 
for <50 = and SO = n, and a minimum kick of 
(352±10)kms- 1 for SO = -tt/2. This will be the subject 
of an upcoming paper by the aut hors. The total recoi l 
velocity also acquires a correction (jSopuerta et al.ll2~006f ) 
for small eccentricities, e, of the form V e = V Icco u (1 + e). 

From lGonzalez et al.l (120061) A = 1.2 x 10 4 kms" 1 and 
B = -0.93. From fits to iHerrmann et all (|2007t) and 
iKoppitz et a il (f2007h we find H = (7 3 ±0.3) x lO^ms' 1 
and from SP2 and iGonzalez et alJ (|2007h (which ap- 
peared after this paper had been submitted) we find 
.ftfcos(0— 0o) = (6, —5.3) xlO 4 kms -1 respectively. Note 
the sign difference showing some of the Vr eC oii dependence 
on 0. 

The in-plane recoil velocity for our SP6 configuration 
is consistent with Eq. |2]) with £ w 88° , however our sim- 
ulation shows strong precession of the spin near merger, 
where a large fraction of the recoil velocity is built up, 
hence it is difficult to accurately determine the spin pa- 
rameters a and to be used in Eq. (|2|). 

4 After completio n of this work we become aware of a new paper 
I Baker et al. ( 2007^) also modeling v± 



A n umber of arguments (|Shapiroll2005t iGammie et al.l 
120041) suggest that spins of supermassive black holes 
(SMBHs) are close to maximal, a/m > 0.8, and per- 
haps as great as 0.99 (jRevnolds et al.ll2005h . Mass ratios 
of binary SMBHs are poorly constrained observation- 
ally, but the luminous galaxies known to harbor SMBHs 
are believed to have experienced a few to several major 
mergers (mass ratios 0.3 < q < 1) over their lifetimes 
(|Haehnelt and Kauffmannl 120021 : lMerrittJl2006f) ; mergers 
with q w 1 were more common in the past. Together 
with the results discussed above, these arguments sug- 
gest that recoil velocities for binary SMBHs in galactic 
nuclei are often of order ~ 10 3 km s~ . Here, we consider 
some consequences of such large recoil velocities. 

Ejection: Central escape velocities from giant ellipti- 
cal galaxies and spiral galaxy bulges are 450 km s _1 < 
v e < 2000 km s _1 , dropping to < 300 km s _1 in 
dwa rf elliptical (dE) and dwarf spheroidal (dSph) galax- 
ies ()Merritt et alJ |2004[ ) . Recoil velocities as large as 
10 3 km s" 1 would easily eject SMBHs from dE and dSph 
galaxies, and in fact there is little evidence for SMBHs 
in these galaxies. However, we note that the mass de- 
pendence of spin-dominated kicks, T^ocoii ~ q 2 , implies 
that recoil velocities might only infrequently be as large 
as in the equal-mass case. If the tight empirical relations 
between SMBH mass and l u minous galaxy properties 
(jFerrarese and Merrittl l2000b iMarconi and Huntl l2003h 
are to be maintained, peak recoil velocities are con- 
strained to be < 500 kms -1 (jLibeskind et al.l 120061) : 
the upper limit on V rcco u is relaxed if most of the 
merger-induced growth of SMBHs took place at low 
redshifts (z < 2) when potential wells were deeper 
(jLibeskind et al.ll2006t ). Ejection of SMBHs from shal- 
low potential wells at high redshift implies a maximum 
z at which the progen itors of present-day SMBHs could 
have started merging (jMerritt et al.ll2004T ). and the exis- 
tence of bright quasars at z « 6 is difficult to reconcile 
with recoil velocities > 10 2 km s _1 unless their SMBHs 
grew very quickly via accretion (Haiman 2004). However 
we note that these and similar conclusions are based on 
an assumed mass ratio dependence for the kicks that is 
invalid if recoil velocities are dominated by spin effects. 

Displacement: The long return time for a SMBH 
ejected near escape velocity implies a substantial prob- 
ability of finding a displaced SMBH in a luminous E 
galaxy, especially if the latter was the site of a recent 
merger dMerritt et al.ll2004HMadau and Quataertl 120041 : 
iVicari et al.ll2006lh A recoiling SMBH carries with it ma- 
terial that was orbiting with velocity v > V^ CCO ii before 
the kick; the size of the region containing this mass is 

772 ~ 1 P C M 8CT 200 (3) 

^recoil \ ^recoil / 

where Mg is the SMBH mass in units of 10 s M® and CT200 
is the nuclear velocity dispersion in units of 200 km s — . 
This radius is sufficient to include the inner accretion disk 
and the broad-line region gas, implying that a kicked BH 
can continue shining for some time as a quasar. Plau- 
sibility of models that explain the "naked" quasar HE 
0450-2958 as an ejected SMBH ()Haehnelt et al.ll2006ft are 
enhanced by the larger kick velocities found here, since 
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v e from the nearby galaxy is > 500 km s~ ; however 
the presence of spectral features associated with narrow 
emission line region gas in the quas ar is still difficult t o 
reconcile with the recoil hypothesis (|Merritt et al.ll2006f) . 

Galaxy cores: The kinetic energy of a displaced SMBH 
is transferred to the stars in a galactic nucleus via dy- 
namical friction, lowering the stellar density and enlarg- 
ing the core before the hole returns to its central lo - 
cation (|Merritt et all 120041 : iBovlan-Kolchin et all 12004ft . 
"Damage" to the co re is maximized for V xeC oi\/v e * 0.7 
(|Merritt et al.ll2004l ). hence the effect could be large even 
in the brightest E galaxies. Observed core masses are 
mostly in the range 0.5 — 1.5M,, consistent with the 
cores having been generated by binary SMBHs without 
the help of kicks ( Merritt! [2006) : however a significant 
fraction have core masses exceeding 2M, suggesting an 
additional contribution from recoils. Anomalously large 
cores in BCG's ("brightest cluster galaxies") might be 
explained in this way, since these galaxies have experi- 
enced the largest number of mergers; this explanation 
would lessen the necessity for altern atives that requir e 
BCGs to contain hypermassive BHs (jLauer et al.l l2006) . 

Jet Directions: We measure both a significant pre- 
merger spin precession of ~ 45° and a change in ex- 
cess of 90° between the initial and final spin vectors, 



verifying the spin - flip phenomenon first discussed by 
iMerritt and Ekersl (|2002l ). Thus, our simulation repre- 
sents a possible model for the merger process responsible 
for generating radio sources with changing jet directions. 
In particular, the highly-spinning large mass black hole 
merging with the smaller mass non-spinning hole is a pos- 
sible model for both the gradual semi-per iodic deviations 
of the jet directions from a straight line (Komossa l2006l ) 
due to precession of the spin (and hence jet) direction, 
and the a brupt change in jet direction forming X-sha ped 
patterns (|Parma et al.1 1 19851 : [Leahy and Parmalll992ft . 
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